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Abstract 

The aspect of self-affine nature of faulting and fracture is widely documented from the data analysis of 
■ both field observations and laboratory experiments. In this direction, Huang and Turcotte have stated 
that the statistics of regional seismicity could be merely a macroscopic reflection of the physical processes 
Q ■ in earthquake source, namely, the activation of a single fault is a reduced self-afhne image of regional 
1) \ seismicity. This work verifies the aforementioned proposal. More precisely we show that the population of: 
(i) the earthquakes that precede of a significant event and occur around its the epicentre, and (ii) the "fracto- 
electromagnctic earthquakes" that are emerged during the fracture of strong entities distributed along the 
activated single fault sustaining the system follow the same statistics, namely, the relative cumulative number 
of earthquakes against magnitude. The analysis is mainly performed by means of a recently introduced 
pH *' nonextensive model for earthquake dynamics which leads to a Gutenberg-Richter type law. We examine the 
O ^ variation of the parameters q and a, which are included in the nonextensive law, for different thresholds of 
magnitude and different radius around the epicentre. Such analysis enhances the physical background of 
the underlying self-afRnity. The parameter q describes the deviation of Tsallis entropy from the extensive 
Boltzmann-Gibbs entropy, and a is the constant of proportionality between the energy released during the 
fracture of a fragment and its size. 

Keywords: Nonextensive Tsallis statistics, Earthquake Dynamics, Preseismic Electromagnetic Emissions, 
self- affinity 



1. Introduction 

The self-affine nature of fracture and faulting has been widely documented from the analysis from both 
laboratory and field observations [1-5]. Characteristically, Huang and Turcotte [2] suggested that the statis- 
tics of regional seismicity could be merely a macroscopic reflection of the physical processes in the earthquake 
(EQ) source. This suggestion, from the perspective of self-affinity, implies that the activation of a single 
fault is a reduced self-affine image of regional seismicity. In this work we attempt to verify the aforemen- 
tioned proposal drawing from a recently introduced model for EQ dynamics [6, 7] which is based on the 
nonextensive Tsallis statistical approach [8-10]. Its theoretical ingredient concerns two rough profiles in- 
teracting via the fragments filling the gap between them. The model leads to a nonextensive Gutenberg & 
Richter (G-R) type formula that describes the frequency distribution of EQs against magnitude including 
two parameters: (i) the entropic index q, which describes the deviation of Tsallis entropy from the standard 
Boltzmann-Gibbs entropy, namely the index q interprets the degree of non-extensivity, that accounts for the 
case of many non-independent, long-range interacting subsystems, and (ii) the physical quantity a, which 
is a constant of proportionality between the energy released during the fracture of a fragment and its size r. 
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The question naturally arises to whether the self-afSne nature of fracture and faulting can be adequately 
explained by the nonextensive G-R type formula. In this direction we examine whether the nonextensive 
G-R type equation fit the following two populations of fracture events: (i) The fractures-EQs that precede a 
significant seismic event occurring in the region which surrounds its epicenter (foreshock activity), (ii) The 
fractures of strong entities which are distributed along the single activated significant fault sustaining the 
system. 

At this point the following question is emerged: "How can we know the sequence of magnitudes of fractures 
of strong entities which are evolved during the activation of a significant fault?". Crack propagation is the 
basic mechanism of material failure. It has been shown that during the mechanical loading, fracture induced 
acoustic and electromagnetic (EM) fields have been observed which allow the real-time monitoring of damage 
evolution in these materials. EM emissions in a wide frequency spectrum ranging from kHz to MHz are 
produced by opening cracks, which can be considered as the so-called precursors of general fracture. The 
radiated EM precursors are detectable both at a laboratory [11-20] and geophysical scale [21-44]. 

An important feature, observed at both scales, is that the MHz radiation precedes the kHz one [29, 31- 
33, 39, 45]. Studies on the small (laboratory) scale reveal that the kHz EM emission is launched in the tail 
of pre-fracture EM emission from 97% up to 100% of the corresponding failure strength [29, and references 
therein]. At the geophysical scale the kHz EM precursors are also emerged in the tail of preseismic EM 
emission, namely, from a few days up to a few hours before the EQ occurrence [32]. Thus, the association 
of MHz, kHz EM precursors with the last stages of EQ generation seems to be justified. 

Based on the above mentioned experimental facts, a two-stage model has been recently proposed con- 
cerning the detected EM emissions in the field [32, 33, 36, 38, 46, 47] suggesting that: 

(i) The MHz EM emission is due to the fracture of the highly heterogeneous system that surrounds 
the fault. More specifically, the MHz EM activity can be attributed to phase transition of second 
order [33], while a Levy walk type mechanism can explain the observed critical state [38]. We note 
that the heterogeneity and long-range correlations are two of the key components that make material 
failure an interesting field for the application of statistical mechanics. The suggestion that ruptures 
of heterogeneous systems is a critical phenomenon has been proposed by numerous of authors, e.g 
[4, 48-53, 53-55] 

(ii) The kHz EM phenomenon is rooted in the final stage of EQ generation, namely, the fracture of backbone 
of entities that prohibit the relative slipping of the two profiles of the fault [32, 33, 36, 37, 39]. Thus, the 
detected precursory sequence of "fracto-EM kHz earthquakes" (see sec. 4.2 ) represents the population 
of fractures that occur during the relative displacement of fault plates. 

We note that the aforementioned proposal has been supported by a multidisciplinary analysis, e.g., in 
terms of extensive and non extensive statistical physics [36-38, 56], information theory, complexity [34, 57], 
laboratory experiments [29, 35, 47], fault modeling [26], fractal electrodynamics [31], self-affinity in fracture 
and faulting [40], nonextensive model for earthquake dynamics [36, 58, 59], and mesomechanics [47]. We note 
that the strong impulsive kHz EM time series shows strong persistent behaviour mirroring a non-equilibrium 
process without any footprint of an equilibrium thermal phase transition. We also note that fracture surfaces 
have been found to be self-affine following the persistent fractional Brownian motion (fBm) model over a 
wide range of length scales, while, the spatial roughness of fracture surfaces H ~ 0.7 has been interpreted as 
a universal indicator of surface fracture, weakly dependent on the nature of the material and on the failure 
mode [60-64, and references therein]. 

In this work, we attempt to provide evidence that the activation of a significant single fault is a reduced 
self-affine image of the regional seismicity, namely, the foreshock neighboring activity associated with the 
activation of the main fault, as follows: We first examine whether the aforementioned nonextensive G-R 
type formula can adequately describe both: (i) the populations of EQs included in different radius around 
the epicenter of a significant seismic event (foreshock activity), and (ii) the populations EM-EQs (see sec. 
4.2) mirroring the fracture of strong entities distributed along the main fault that sustain the system. 
Furthermore, we focus on the variation of parameters q and a included in the nonextensive formula for both 
populations using different thresholds of magnitudes. In this direction, two well documented cases that fulfil 
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the needs of such analysis, with available data for both the seismicity and preseismic kHz EM emissions 
observed prior to large EQs, are: (i) the case of Athens EQ (M ~ 5.9) occurred on 07-Sep-1999, in Greece 
and (ii) the case of L'Aquila EQ {M = 6.3) occurred on 06-Apr-2009 in Italy. 

This paper is organized as follows: In Sec. 2, we describe the basic principles of Tsallis nonextensive 
statistical mechanics and the nonextensive model for EQ dynamics. In Sees. 3 and 4, we apply nonextensive 
analysis on both the seismicity and preseismic kHz EM emissions, for the case of L'Aquila EQ and Athens 
EQ, respectively. In Sec. 5 we analyze our results providing further arguments that support the self-affine 
nature of fracture and faulting. Finally in Sec. 6, we summarize the key findings. 

2. Theoretical background 

The aim of statistical mechanics is to establish a direct link between the mechanical laws and classical 
thermodynamics. Within that context, "extensivity" is one of the crucial properties of the Boltzmann-Gibbs 
entropy (Sb-g) expressing the proportionality with the number of elements of the system. Sb-g satisfies 
this rule, if the subsystems are statistically (quasi-) independent, or typically if the correlations within the 
system are essentially local. In such cases the system is called extensive. However, in the cases where the 
correlations may be far from negligible at all scales, the Sb-g is "nonextensive". Inspired by multi-fractal 
concepts, Tsallis [8] proposed a generalization of the B-G statistical mechanics, by introducing an entropic 
expression characterized by an index q which leads to a nonextensive statistics: 

where pi are the probabilities associated with the microscopic configurations, W is their total number, q is 
a real number, and k is Boltzmann's constant, q — 5- 1 corresponds to the standard extensive B-G statistics. 
Indeed, using pj* = e^''"^-' '"^^'^ ~ 1 + (q — l)ln(pi) in the limit q — > 1, we obtain the standard B-G 
entropy: 

w 

Si ^ -kY^p.Hp.) (2) 

4=1 

The entropic index q characterizes the degree of non-additivity refiected in the following pseudo-additivity 
rule: + B) ^ Sg{A) + Sg{B) + (1 - q)SgiA)Sg{B). 

The cases q > 1 and q < I, correspond to sub-additivity, or super-additivity, respectively. The parameter 
q itself is not a measure of the complexity of a time series. It measures the degree of nonextensivity of the 
corresponding system. A metric of the dynamic changes of the complexity of a system is the time variations 
of the Tsallis entropy for a given q (Sq). 

2.1. A fragment- asperity model for earthquakes coming from a nonextensive Tsallis formulation 

The best known scaling relation for EQs is the Gutenberg & Richter (G-R) magnitude-frequency rela- 
tionship [65], given by: 

\ogN{> m) ^ a -hm, (3) 

where iV(> m) is the cumulative number of EQs with a magnitude greater than m occurring in a specified 
area and time. Parameters h and a are constants. 

A nonextensive model for EQ dynamics consisting of two rough profiles interacting via fragments filling 
the gap has been recently introduced by Sotolongo-Costa and Posadas (SCP) [6]. More recently, the afore- 
mentioned model was revised by Silva et al. [7], where two crucial ingredients were employed. They use 
revised definition for the mean values in the context of Tsallis nonextensive statistics that was achieved in 
the study of Abe and Bagci [66]. Moreover, Silva et al. proposed a new scaling law, e o< r'^ , between the 
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released relative energy e and the size r of fragments. Finally, their approach leads to the following G-R 
type law for the magnitude distribution of EQs: 



G(> M) 
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where, N is the total number of EQs, N{> M) the number of EQs with magnitude larger than M, and 
M w log(e). The parameter a is the constant of proportionality between the energy, e, and the size of 
fragment, r. This is not a trivial result, and incorporates the characteristics of noncxtensivity into the 
relative cumulative number of EQs against magnitude. The entropic index q describes the deviation of 
Tsallis entropy from the traditional B — G one, expressing the long-range correlations developed in the 
system. 

The g-parametcr included in the non-extensive formula (Eq. (4)) is associated with the b parameter of 
Gutenberg & Richter formula (Eq. (3)), by the relation [67]: 



6 = 2 X 



2-q 



(5) 



Note that the nonextensive model of Eq. (4) has been successfully applied for several different regions 
under study incorporating the characteristics of nonextensivity statistics into the detected EQs [7, 68-72] . 
Values of q ^ 1.6 — 1.8 [7, 68-70, 73] seem to be universal, in the sense that different data sets from different 
regions of the globe indicate a value for the nonextensive parameter lying in this interval. In the following 
sections the two aforementioned cases of strong EQs are examined, for both the seismicity and the observed 
preseismic kHz EM emissions. 



3. The case of L'Aquila 2009 earthquake 

We recall that the prospect in this work is to examine whether the self-affine nature of fracture and 
faulting can be adequately explained by the nonextensive Eq. (4). Eq. (4) is directly connected to the 
traditional G-R law (Eq. 3), above some magnitude threshold through Eq. (5), which in turn leads to 
the power-low distribution of magnitudes, expressing the fractal nature of the system under study [74-78]. 
Thus in the prospect to verify the aforementioned suggestion, in the following sections we examine both the 
seismicity (foreshock activity) and the observed preseismic kHz EM emissions related to large EQs, in the 
prospect to verify the aforementioned suggestion. 

3.1. Analysis of seismicity (L'Aquila case) 

The analysis here is focused on the very shallow strong (Mw = 6.3) earthquake occurred on 6 April 2009 
in L'Aquila, Central Italy at 01:32:39 UTC. Characteristically, Papadopoulos et al. [79], reported that from 
the beginning of 2006 up to the end of October 2008 no particular earthquake activation was noted in that 
seismogenic area. On the contrary, from 28-October-2008 up to 27-March-2009 the seismicity was in the state 
of weak foreshock activity, and dramatically increased 10 days prior the main event (from « 26-Mar-2009) 
[79]. Drawing from these evidence, the period from 28-Oct-2008 00:00:00 up to 6-April-2009 01:32:00 was 
selected as the most appropriate time- frame in order to examine the sclf-afhne behaviour of the EQs included 
within different radii around the EQ epicenter. The Italian EQ catalogue which is available on the website of 
the Istituto Nazionale di Geofisica e Vulcanologia (INGV: http://bollettinosismico.rm.ingv.it), was 
used. 

Herein, we first examine whether the aforementioned nonextensive G-R type formula (Eq. (4)) can 
adequately describe the populations of EQs included within different radii around the epicenter of L Aquila 
EQ. In this direction Eq. (4) was used to fit the seismic data of six different selected geographic areas 
around the EQ epicenter, for the period under study: 0-400km, 0-300km, 0-200km, 0-100, 0-50km and 
0-30km respectively, as shown in Fig. 1. The Levenberg-Marquardt (LM) fitting method [80, 81] was 
used for optimal fitting. From Fig. 1, we observe that Eq. (4) provides an excellent fit on the relative 
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cumulative number of magnitudes {G{> M)) related to the populations of EQs included in each one of the 
six selected geographic areas. The noncxtensive g-parameter varies between [1.644 1.694] with a relative 
small standard error ranging between [0.002 0.003]. 



Silva et. al. formula Silva et. al. formula Silva et. al. formula 
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(d) 0-lOOkm (e) 0-50km (f) 0-30km 

Figure 1: We use Eq. (4) to fit the foreshock seismic data in terms of the relative cumulative number of 
EQs included in six different geographic areas around the L'Aquila EQ epicenter. 

Secondly, in the prospect to examine the behaviour of parameters q and a included in the nonextensive 
formula, different thresholds of magnitudes cut-offs (Mc) were applied using an increasing step of 0.1. For 
each step, formula (4) along with the use of Levenberg-Marquardt (LM) fitting method, were used to fit 
the seismic data in terms of the relative cumulative number of EQs. A minimum number of 20 events was 
considered as a criterion for the statistical completeness of each calculation. The derived parameters q and 
a, were graphically placed on to a common x-axis chart as shown in Fig. 2. It is observed that both the 
nonextensive g-parameter and the energy a, present similar behaviour, for all the selected geographic areas 
around the EQ epicenter. More specifically, the nonextensive parameter q (depicted with black bullets) 
remains relative stable with minor increment at intermediate thresholds of magnitudes (Mc 1.5 — 2.5). 
The characteristic value of the volumetric energy density a (depicted with blue- rhombuses ) , increases at 
higher magnitude thresholds (Mc). 
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Figure 2: (a) Variation of noncxtcnsivc parameter q (see black-bullets) and the volumetric energy density a 
(see blue- rhombuses) , for different thresholds of magnitudes of the detected EQs included in the period from 
28-Oct-2008 00:00:00 up to 6-April-2009 01:32:00, for different geographic areas around the EQ epicenter: 
(a) 0-400km (b) 300km, (e) 0-200km (d) 0-lOOkm (e) 0-50 and (f) 0-30km around the epicenter of L'Aquila 
EQ. 



As we have seen so far, the analysis of seismicity at different geographic scales revealed similar footprints 
for both the parameters included in the nonextensive formula. The same behaviour is expected to be found 
at lower scales of fracturing, namely, the preseismic kHz EM anomalies observed prior to L'Aquila EQ. 
Thus, in the following subsection, we examine these kHz EM emissions in terms of the nonextensive model. 

3.2. Analysis of preseismic kHz EM emissions observed prior to L'Aquila EQ 

This subsection focuses on the kHz EM emissions detected from a measurement station which has been 
installed and operating at a mountainous site of Zante Island in the Ionian Sea (Western Greece). During 
the aforementioned 10-day foreshock activity, well documented [39, 45], preseismic kHz EM anomalies were 
observed on 4-Apr-2009, two days after the MHz EM anomalies of 2*'' Apr, verifying the two-stage model 
described in Sec. 1. More precisely, the detected anomalies followed the temporal scheme listed below: 

(i) The MHz EM anomalies were detected on 26 March 2009 and 2 April 2009. 

(ii) The kHz EM anomalies emerged on 4 April 2009. 

(iii) The ULF EM anomaly was continuously recorded from 29 March 2009 up to 2 April 2009. 

In Fig. 3, the top charts of each sub-figure refer to the three observed magnetic field strengths, recorded 
from the 10 kHz (NS,EW,V) sensors, on 4-Apr-2009. Note that the EQ occurred on 06-Apr-2009 01:32:39 
UTC as shown by the black-arrow. The selected red parts refer to the period where the kHz EM anomalies 
observed, which in turn have been well justified for their seismogenic origin from recent studies [39, 45]. 

The notion of electromagnetic earthquake: We regard as amplitude ^ of a candidate "fracto-electromagnetic 
fiuctuation" the difference Afem(ti) = ^{tj) ~ ^noise, where Anoise is the background (noise) level of the 
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EM time scries. We consider that a sequence of k successively emerged "fracto-electromagnetic fluctua- 
tions" Afem{ti), i = 1, . . . , fc represents the EM energy released, e, during the damage of a fragment. We 
shall refer to this as an "electromagnetic earthquake" (EM-EQ). Since the squared amplitude of the fracto- 
electromagnetic emissions is proportional to their power, the magnitude m of the candidate EM-EQ is given 
by the relation: 

m = loge - log [A/e,„(<j)]^^ (6) 

Herein, as in the case of seismicity, we first examine whether the aforementioned nonextensive G-R type 
formula (Eq. (4)) can adequately describe the populations of EM-EQs included in the preseismic kHz EM 
time series. As shown in the middle charts of Fig. 3, we used Eq. (4) to fit the EM data in terms of 
the relative cumulative number of electromagnetic earthquakes G{> M), included in the period that refers 
to the red part of these signals. It is observed that Eq. (4) provides a satisfactory fit to the preseismic 
kHz EM experimental data associated with the L'Aquila EQ for all the three recorded components. The 
nonextensive g-parameter varies between [1.818 ^ 1.838] with a relative small standard error ranging between 
[0.001 - 0.002]. 




Figure 3: Three observed magnetic fields recorded on 02-Apr-2009 16:00:00 to 06-Apr-2009 01:32:39. The 
selected red parts of each recorded signal refer to the period where the kHz EM anomalies observed and 
have been well justified for their seismogenic origin ([32, 33, 36, 38]). The middle charts of each subfigure, 
refer to the fitting of the EM data in terms of the relative cumulative number of EM-EQs G(> M) (see 
Eq. (4)), included in the period under study (red-part). The bottom charts of each sub-figure, refer to the 
corresponding variation of the nonextensive g-parameter (black-bullets) and the volumetric energy density 
a (blue-rohmbs), for different thresholds of magnitude cut-off (Mc). 

Secondly, we apply Eq. 4, on the relative cumulative number of EM-EQ contained in the preseismic kHz 
EM emissions in the prospect to examine whether the variation of nonextensive g-parameter the volumetric 
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density a, also present the same behaviour with that of seismicity (see See. 3.1). Different thresholds of 
magnitude cut-offs (Mc) were applied on the EM-EQs contained in eaeh one of the aforementioned channels 
using an increasing step of 0.1. For each step, Eq. (4) along with the use of LM method, was used to 
fit the EM data in terms of the relative cumulative number of EM-EQs contained in each channel. In 
Fig. 3, the bottom charts of each sub-figure, depict the variation of the nonextensive g-parameter (black 
bullets) and the volumetric energy density a (blue rhombuses), for different thresholds of magnitudes (Mc). 
It is observed that the variation of the nonextensive g-parameter remains relative constant with minor 
decrement at higher thresholds of magnitudes. On the contrary, the energy density a mirrors the behaviour 
of g-parameter, presenting a relative increment at higher thresholds. 

These latter results evidently show the similarity with the results obtained in Sec. 3.1, indicating the 
self-affine nature of fracture and faulting, from the large scale of foreshock seismicity, to the fault-generation 
of a single EQ in terms of preseismic kHz EM emissions. However, in order to ensure the consistency of the 
results, in the following sections, we further examine the case of Athens 1999 EQ. 
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4. Additional evidence of the self-afRnity of EQs: The case of Athens 1999 earthquake 

4-.1. Analysis of seismicity (Athens case) 

Wc focus on the case of Athens EQ {Mw = 5.9) occurred on 07-Scp-1999 11:56:50 UTC, namely, the 
period from 17-Aug-1999 00:01:39.80 up to 07-Sep-1999 01:56:49. The Greek EQ catalog was used, as 
provided by the website of the Institute of Geodynamics of the National Observatory of Athens (www.gein. 
noa.gr). Note that Athens EQ occurred very shortly after the major 17/8/1999, Mw = 7.4, EQ which took 
place on Izmit, Turkey, approximately 650 km North-East of Athens. It has been shown that the Izmit EQ, 
was followed immediately by small earthquakes occurred throughout much of continental Greece [82]. This 
explains the selected period. 

Four different geographic areas were selected around the Athens EQ epicenter for the period under 
study: — 400fcm, — SOOfcm, — 200^™ and — 160fcTO respectively. In Fig. 4, we present the fitting 
of Eq. (4) applied on the relative cumulative number of magnitudes (G(> M)) for the populations of EQs 
included in each selected geographic area. It is observed that the nonextensive q-parameter varies between 
[1.532 ^ 1.664] with a relative small standard error ranging between [0.003 ~ 0.004]. 



Silva et. al. formula Silva et. al. formula 
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Figure 4: We use Eq. (4) to fit the seismic data in terms of the relative cumulative number of EQs included 
in four different geographic areas around the Athens EQ epicenter. 



Applying the same method as that described in Sec. 3, in Fig. 5, we present the variation of nonextensive 
parameter q (see black-bullets) and the volumetric energy density a (see blue rhombuses), using different 
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thresholds of magnitudes. It is observed that both the nonextcnsive g-parameter and the energy a, present 
similar behaviour, for all the selected geographic areas around the EQ epicenter. More specifically, the 
nonextcnsive parameter q (depicted with black bullets) remains relative stable with minor decrement at 
higher thresholds of magnitudes, while the characteristic value of the volumetric energy density a (depicted 
with blue rhombuses), increases at higher magnitude thresholds (A'/c)- 
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Figure 5: Variation of nonextcnsive parameter q (see black-bullets) and the volumetric energy density a 
(see blue- rhombuses) , for different thresholds of magnitudes of the detected EQs included in the period from 
17-Aug-1999 00:01:39.80 up to 07-Sep-1999 01:56:49, for the range (a) 0-400km (b) 0-300km, (c) 0-200km 
and (d) 0-160km around the Athens EQ epicenter. 



4-2. Analysis of preseismic kHz EM emissions observed prior to Athens EQ 

The analysis here, is focused on the well documented [32, 33, 36, 38] kHz EM activity (Fig. 6, upper 
sub-charts) observed before the Athens EQ. The the six observed magnetic field strengths recorded by the 3 
kHz (NS, EW and V) and 10 kHz (NS, EW and V) sensors, are analyzed, in terms of the nonextcnsive Eq. 
(4). The top charts of each sub-figure depicted in Fig. 6, refer to the six observed magnetic fields recorded 
on the period from 28-Aug-1999 00:00:00 to 08-Sep-1999 00:00:00. Note that the Athens EQ occurred on 
07-Sep-1999 07-Scp-1999 01:56:50 UTC as shown by the black-arrow. The selected red parts refer to the 
period where the kHz EM anomalies observed, which in turn have been well justified for their seismogenic 
origin from recent studies [32, 33, 36, 38]. 

We note that in a recent study [59] the six observed magnetic fields have been analyzed in the context 
of nonextcnsive Tsallis statistics. The analysis showed that the g-parameter lies within the interval q £ 
[1.78 — 1.82] with a relative small standard error ranging between [0.0003 ^ 0.0006], for all the recorded 
channels. Herein, applying the same method as that described in Sec. 3.2, we further investigate the 
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behaviour the noncxtensive g-parameter and the volumetric energy density a from the perspective of self- 
affinity. More precisely, the bottom charts of each sub-figure contained in Fig. 6, depict the variation of the 
nonextensive g-parameter (black bullets) and the volumetric energy density a (blue rhombuses), for different 
thresholds of magnitudes (Mc). 

From Fig. 6, it is observed that the variation of the nonextensive (/-parameter remains relative constant 
with minor decrement at higher thresholds of magnitudes. On the contrary, the energy density a mirrors 
the behaviour of g-parameter, presenting a relative increment at higher thresholds. These results evidently 
show the similarity with the results obtained in Sec. 4.1, indicating the self-affine nature of fracture and 
faulting, from large-scale seismicity, to the fault-generation of a single EQ in terms of preseismic kHz EM 
emissions. In the next section we further provide further analysis and arguments that support the findings 
of this work. 



28-Aug-1 999 00:00:00 to 08-Sep-1 999 00:00:00 




(d) 10 kHz NS (c) 10 kHz EW (f) 10 kHz V 



Figure 6: Six observed magnetic fields recorded on the period from 28-Aug-1999 00:00:00 to 08-Sep-1999 
00:00:00. The selected red parts of each recorded signal refer to the period where the well justified for 
their seismogenic origin kHz EM anomalies have been observed ([32, 33, 36, 38]). The bottom charts of 
each sub-figure, refer to the corresponding variation of the nonextensive g-parameter (black-bullets) and the 
volumetric energy density a (blue-rhombuses), for different thresholds of magnitude cut-off (Mc)- 
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5. Analysis of results - Discussion 

Herein we recall that the main target of this work is to verify the self-afhne nature of fracture and faulting 
in the framework of nonextensive Tsallis statistics; we examined whether it can be adequately explained by 
the nonextensive G-R type formula (Eq. 4). In this direction we note that the traditional G-R formula leads 
to the power- low distribution of magnitudes expressing the fractal nature of the system under study [74-78] . 
On the other hand, Eq. (4) is directly connected to the traditional G-R law (Eq. 3), above some magnitude 
threshold through Eq. (5) [58, 59, 67, 71, 83]. Moreover, since the released energy e is proportional to 
fragment size (e oc r^) [7], it is reasonable to assume that the magnitudes of EQs which are rooted in 
the fracture of the population of fragments-asperities filling the space between fault planes, also follow a 
power-law distribution. 

The formula was mainly applied on both: 

(i) the populations of EQs included in different radius around the epicenter of a significant seismic event 
(foreshock activity), and 

(ii) the populations of EM-EQs which mirror the fracture of strong entities distributed along the main 
fault that sustain the system. 

Indeed, analysis on both the populations of EQs included in foreshock activity and the populations of 
EM-EQs included in preseismic kHz EM emissions (see Figs. 1,3 and 4) has provided similar footprints 
of nonextcnsivity verifying the aforementioned relation. More specifically, for the case of L'Aquila EQ the 
nonextensive g-parameter ranges between [1.644 ^ 1.694] for the seismicity and between [1.818 ~ 1.838] 
for the preseismic kHz EM emissions. For the case of Athens EQ, the nonextensive g-paramcter ranges 
between [1.532 ^ 1.664] for the seismicity and between [1.78 ^ 1.82] for the preseismic kHz EM emissions. 
These calculated nonextensive g-parameters are in general agreement with those obtained from several 
independent studies, related to seismicities generated in various large geographic areas, involving the Tsallis 
nonextensive framework [7, 68-70, 72]. We indicatively mention the study of Telesca [70], who found a value 
of {q = 1.742) by examining the preseismic activity of L Aquila EQ, included in the period from 30-Mar-2009 
to 06-Apr-2009. 

Focusing on the method for the estimation of the behaviour of parameters q and a, by using different 
magnitude thresholds (Mc), similar behaviour has also been found for both the populations of EQs included 
in foreshock activity and the populations of EM-EQs included in preseismic kHz EM emissions. Charac- 
teristically, from Figs. 5, it is observed that the nonextensive parameter q and energy a, present similar 
behaviour, for all the selected areas around the Athens EQ epicenter. More specifically, it is observed that 
the nonextensive parameter q (depicted with black bullets) remains relative stable for different magnitude 
thresholds. The situation changes for even larger thresholds of magnitude, where a relative decrement is 
observed. Such behaviour is also observed from the analysis of preseismic kHz EM emissions as shown in Figs 
3 and 6. This prospective decrement can be explained by the fact that the larger the magnitude threshold 
the larger the number of the omitted EQs. The absence of the small fractures along with the corresponding 
redistribution of stresses, contributes to the decrement of the correlation length during the fracture process 
[84]. The smaller magnitude threshold is the one that governs the overall system (e.q for Mc = 0). It should 
also be noted that although the nonextensive parameter q decreases at higher magnitude thresholds it still 
remains high veriiying the strong correlations that have been developed. Concerning the volumetric energy 
density a (depicted with blue-rhombuses), its characteristic value, increases at higher thresholds of magni- 
tude {Mc). Note that according to the fragment-asperity model (SCP), a is the coefficient of proportionality 
between fragment size and released energy [6, 7]. This evidence is consistent with the hypothesis that larger 
EQs are rooted in larger and stronger entities [36, 58] that sustain the system. 

A characteristic differentiation that should also be discussed here, mainly observed in the case of L Aquila 
analysis of seismicity (see Fig. 2), is the relative inclement of q parameter at higher magnitude thresholds 
as opposed to the case of Athens EQ, where a relative decrement is observed (see Fig. 5). This is also not 
a trivial result because if we get a closer look to Athens case in Fig. 5, this variation is mainly observed for 
larger thresholds of magnitude in contrast to the case of L Aquila EQ where the increment is observed at 
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intermediate thresholds of magnitude. This evidence further enhances the suggestion that the smah fractures 
along with the corresponding redistribution of stresses, contributes to the increment of the correlation length 
during the fracture process [84] . Indeed we calculated the mean magnitude of the EQs included in the area 
of 0-400km around the epicenter, for both the examined periods and has been found to be: M = 1.60 for 
the case of Italy EQ and M = 3.00 for the case of Greece. Furthermore, Fig 7 shows the distribution of EQ 
magnitudes that correspond to the areas under study. It is clearly observed that in the case of the Italian 
territory, EQs with lower magnitudes are more frequent in contrast to the Greek territory. 
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Figure 7: Distribution of EQ magnitudes included in the area of 0-400km around the Italian EQ (red bars) 
and Greek EQ (blue bars), correspondingly. 

Along with these results, another feature that is also observed as approaching to the EQ epicenter, is 
that the smaller the area under study around the epicenter, the higher the nonextensive g-parameter. This 
evidence further verifies the increased correlations developed as approaching to the EQ epicenter. On the 
contrary the energy a, presents a relative decrement that mirrors the behaviour of g-parameter. However, 
this mirroring behaviour between q and a, is still an open issue for the scientific community and also out of 
the scope of this study. 

6. Conclusions 

Building on the perspective of the self-affine nature of fracture and faulting theory, this work supports 
the hypothesis that the statistics of regional seismicity is a macroscopic reflection of the physical processes 
in the earthquake source. This suggestion implies that the activation of a single fault is a reduced self-affine 
image of regional seismicity. We used a recently introduced nonextensive Gutenberg & Richter type formula 
which describes the EQ dynamics and includes two parameters: the entropic index q, which describes the 
deviation of Tsallis entropy from the standard Boltzmann-Gibbs entropy, and the physical quantity a, 
which characterizes the energy density. Focusing on two cases of large EQs, we found similar behaviour 
of the corresponding variation of the parameters q and a, which are included in the nonextensive law for 
different thresholds of magnitudes applied on: (i) the population of EQs contained in different radii around 
the epicentre of a preceding large EQ and (ii) the population of " fracto-clectromagnetic earthquakes" that 
are emerged during the fracture of strong entities distributed along the activated single fault sustaining 
the system. Analysis revealed that these two populations follow the same statistics, namely, the relative 
cumulative number of EQs against magnitude. Furthermore the results are further supported from recent 
studies in terms of the traditional Gutenberg-Richter law [74-78]. Such analysis enhances the physical 
background of the underlying self-affinity. 
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